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Abstract 

The frequency-dependent exchange-correlation potential, which appears in 
the usual Kohn-Sham formulation of a time-dependent linear response prob- 
lem, is a strongly nonlocal functional of the density, so that a consistent 
local density approximation generally does not exist. This problem can be 
avoided by choosing the current- density as the basic variable in a general- 
ized Kohn-Sham theory. This theory admits a local approximation which, for 
fixed frequency, is exact in the limit of slowly varying densities and perturbing 
potentials. 
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Gross and Kohn (GK) JT| have applied the time dependent density functional theory 

(TDFT) of Runge and Gross, to the case of the linear response of interacting electrons 

in their ground-state to a time-dependent potential Vi(r, u)e~ llJt . The objective was the 

determination of the linear density response, ni(r, uj)e~ tu)t . They reduced the problem to a 

set of self-consistent single particle equations, analogous to the Kohn-Sham equations for 

time-independent systems |y, with an effective potential of the form 

e.ff/-> \ \ fnAr'^uj) , , . ... 

v x JJ (r,u) = vi(r,u) + / — —-dr + v lxc (r,u); (1) 

J \r — r\ 

the exchange-correlation (xc) potential Vi xc (r, u) is linear in ni(r,u), 

vi xc (r, (*})=! fxc(r, r"; u)^^, ^dr 1 , (2) 

and the kernel f xc (r, f ; u) is a functional of the unperturbed ground state density no(r). 

In the spirit of the local density approximation (LDA) for static and quasi-static problems 
0] , they then considered the case where both n and n\ are sufficiently slowly varying 
functions of r. As f xc is of short range for a homogeneous system, they proposed the 
following plausible approximation for systems of slowly varying no{f): 

f xc (r,f>;uj) ~ f h xc (\r-7\,u-n {r)). (3) 

The superscript h refers to a homogeneous electron gas and the function f£ c is a property 
of the homogeneous electron gas [|].|5| . 

However it was noted later by Dobson that the approximation (|3|), when applied to 
an electron gas in a static harmonic potential |&x 2 and subjected to a uniform electric field, 
V\{f,u) = —E ■ fe~ luJt , violates the so called harmonic potential theorem (HPT), related 
to the generalized Kohn's theorem 0, according to which the density follows rigidly the 
classical motion of the center of mass: 7ii(r, uS) = Vn (r) • R C m{u). This raised serious 
questions about the validity of the approximation (|3|). Dobson observed that one could 
satisfy the HPT by requiring that the GK approximation (^[) be applied in a frame of 
reference moving with the local velocity of the electron fluid. The xc potential obtained by 
this construction is a functional of the current- density as well as the density ||. 
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Further light on the problem with approximation (|3|) was thrown by Vignale's observa- 
tion || that the covariance of the time-dependent Schrodinger equation under transformation 
to an accelerated frame of reference requires the total force exerted on the system by the 
exchange- correlation and Hartree potentials to vanish, in agreement with the third law of 
Newtonian mechanics. This implies that the exact f xc must satisfy the sum rule 



J f xc (f, 7; u)V'no^W = V*w(r-), (4) 
where VQ XC (r) is the static xc potential. This sum rule is violated by eq. (H). More generally, 



one can deduce |T0[ that f xc (r, f , to) for a non-uniform system is of long range in space and a 
nonlocal functional of the density distribution. These results indicate that, contrary to more 
optimistic expectations, a local-density approximation for time- dependent linear response 
in general does not exist as long as one insists on describing dynamical exchange-correlation 
effects in terms of the density. 

In this paper, we want to demonstrate however that a local approximation for the time- 
dependent linear response theory can be constructed in terms of the current density. We 
consider the linear current response ji(r,uj)e~ luJt of interacting electrons in their ground- 
state to a time- dependent vector potential ai(r, u)e~ luJt . This problem includes, as a special 
case, the scalar potential problem studied by GK, because any scalar potential v\{f,uS) can 
be gauge-transformed to a longitudinal vector potential ai(r, uS) = Vvi(f,u)/iu, and the 
density response can be calculated from the current response using the continuity equation 
ni(r, uj) = V-ji(r, uj)/iuj. As usual, we express the exact induced current as the response of a 
non-interacting reference system (the "Kohn-Sham" system) to an effective vector potential 
a 1 — a\ + am + ai xc : 

j u (r, u) = Jj2 XKS,ij(r, f, u) ■ '(f ', uj)d? ', (5) 
j 

where 



a 1H {f, u) = j— j / — ^-(V • ji(r , u))dr (6) 



is the longitudinal vector potential corresponding to the dynamic Hartree potential of eq. (HI), 
and 



J fxcir, • ji(f , ,u)df' (7) 



is a linear functional of the current. The 3x3 tensor kernel f xc is defined as 

f xc ,ij (r , = x~Ks,ij (r,r,uj)- xlj if, r , u) - V» ^-^ Vj , (8) 
where x an d Xks are the current response tensors [|l^] of the interacting system and the 



"Kohn-Sham" noninteracting system respectively. The latter is given as 

XKS,ij{r,r ,u) = 5{r - r )8 ij + — 2^(/ Q - /g) ? x , • ' 9 

where ip a {r) are the solutions of the static Kohn-Sham equation, with eigenvalues e a . (We 
have put e = c = 1). For a homogeneous electron gas of density function of r — f, 

which can be Fourier-transformed to f£ c (A;,a;). In the limit of small wave vector k , which, 
at fixed frequency u, means k « kp,uj/vp (kp = Fermi momentum, vp = Fermi velocity) 
|T3T1 , f xc has the form 

= ^[fLMmkj + ft T (u,n)(k% - kikj)], (10) 

where the factor 1/u 2 has been put in evidence so that the function f^ cL (uJ, n) coincide with 
the k — > limit of f xc (k, u, n) introduced by GK (see eq. @). Both f^ cL {oj, n) and /^(o;, n) 
can, in principle, be computed from the current response function of the homogeneous 
electron gas (see eq. (^)). While some aspects of these functions are known P|,|P|, further 
work is needed for a complete evaluation. 

From eq. flTOp we can obtain the xc potential in real space for a homogeneous electron 
gas subjected to a perturbation which is slowly varying on the scales kp 1 (~ interelectron 
distance) and Vp/uj (~ distance travelled by an electron during a period of the perturbing 
field): 

%J?>o>) = -4{^ n)V • /i(r»] - V x [ft T {u,n)V x Uf,uj)]}. (11) 



Next, in the spirit of the local density approximation, let us consider a system in which 
the unperturbed static density n (r) is slowly varying on the scales of the local kp 1 and 
vf/uj, but not necessarily on the scale of k~ l , the wavelength of the perturbing field. 

The simplest case is that of a periodically modulated electron gas, where the unperturbed 
density, given by 

noif) = n(l + 27 cos (q ■ r)), (12) 

is not only slowly varying (i.e., q << kF,u/vp), but, also, almost uniform, i.e. 7 << 1. 

We shall compute the exact f xc for this system, to first order in 7. It is convenient to 
represent f xc in momentum space. Translational invariance of the uniform density n gives, 
to first order in 7, f xc (k, k, uj) = f xc (k, uj) (see eq. fliTf)) and f xc (k + mq, k, uj) = for integers 
m with \m\ > 1. It remains to calculate the matrix element i xc {k + q, k,u). 

This is facilitated by two exact identities which follow from the transformation of the 
time dependent Schrodinger equation in accelerated reference frames |[L5|| . The physical 
content of these identities is that the total force and the total torque exerted by the Hartree 
and xc potentials on the system must vanish, in accordance with Newton's third law. In the 
specific instance considered here, the identities take the form 

lim k ^ f XC:ij (k + q, k, uj) = ~^[5f xcL qiqj + flM^ij ~ QiQj)), (13) 

and 

lim k ^ eijk-^-f xc ,ij(k + q, k, uj) = ^[Sf xcL - 3f xcT ] ^Qk, (14) 

j,k OK k 10 k 

where is the Levi-Civita tensor and we have suppressed uj and n in f^ cL and f xcT ; and 

5 f n) = f* cL (uj, n) - f h xcL {uj = 0, n). (15) 

These identities are valid for k, q << kp,uj/vF and to first order in 7. The condition on q 
has been used to replace f xcL (q, uj) and f xcT (q, uj) by their small q limit PB| . A third limiting 
form of f xc ,ij is obtained from the usual Ward identity 



lim q ^ i xc (k + q, k, u) = ^ n ^ xc ^^ ,n \ (16) 

Finally, we note the symmetry relation 

%tc,%j{k + q, k, uj) = f xc ,ji{-k, —k -q,uj). (17) 

Equation ( |TTD requires that, for small k and q, f XCj ij(k + q, k,oj) be expressible as a linear 
combination of q^j, q 2 5 i j 1 (ki + qi)kj, k{(kj + qj), and k ■ (k + q)Stj. The equations [TB|, [TJ], and 



|To] , completely determine the coefficients of this combination. Thus, to first order in 7, we 
obtain 

W£ + q, k, 00) = -^ 2 {{5f h xcL - f h xcT ) mj + f* cT q% - n^k ■ (k + q)5 t , 

+A(n,uj)(ki + qi)kj - B(n,u)ki(kj + qj)}, (18) 

where A(n,u>) = [n(2df^ cT /dn - dftJdn) + 3f* cT - 8ft L ], and B{n,u) = [ndf£ cT /dn + 
3fxcT ~ &fxcL\- The essential point is that eq. (0) is analytic for small k and q. It is 
this feature that will enable us to construct a local approximation for a\ xc in terms of the 
current density (By contrast, the off-diagonal component f xc {k + q, k,u) in the usual time 
dependent density functional theory, has a singularity of the form k ■ q/k 2 when k — > at 
finite q. This small-fc singularity reflects the already mentioned long-rangedness of f xc and 
the consequent nonlocality of the description of xc effects in terms of the density.) 

Let us now return to the more general problem of determining a\ xc for a system whose 
density is slowly varying in the sense |Vno(r)|/no(r) << kp(r) and uj/vjp(r), but may have 
large global deviations from uniformity. We assume that a lxc of such a system can be 
expanded in a power series of gradients of the local unperturbed density. The most general 
form of Si xc containing up to two gradient operators is then a linear combination of the 
following terms: J x , Vn x j 1; V x ji, V(V ■ ji), V x (V x Ji), V(Vn Q • ji), V x (Vn x 
ji), Vn (V- ji), Vn x(Vxji), (Vn - V)ji, Vn (Vn -ji), Vn x (Vn x j x ) with coefficients 
that are functions of n (r) and uj. This gradient expansion is applicable, in particular, to 
the periodically modulated electron gas of eq. (|12|) , in which case it must yield the same 
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perturbative (in 7) result that one obtains from the combination of eqs. [7|, [10], and [18|. Thus, 
by requiring agreement between the gradient and the perturbative expansions, we are able 
to determine the coefficients of all the terms appearing in the gradient expansion, except 
the last two, which are of order |Vno| 2 , or 7 2 in the periodically modulated electron gas. 
We note that the part of a lxc which is determined by this procedure satisfies the HPT up to 
corrections of order |Vno| 2 . To fix the coefficients of the second order terms, an additional 
condition must be imposed. We require that, if the system is subjected to an external 
field, which causes it to translate rigidly as a whole, i.e., jx(r,u) = no(f)v((j), then a\ xc 
must reduce to — V[/ xc z,(c*; = 0)V • ji(r, uj)]/u 2 , a rigid displacement of the static exchange- 
correlation potential ||. In the special case of electrons confined by a harmonic potential 
and subjected to a uniform electric field, this condition guarantees that the HPT is satisfied. 

The final result of our analysis, obtained after lengthy but elementary manipulations, 
is the complete form of the local approximation for ai xc (f, uf) up to second order in the 
gradient expansion: 

ax xc (r, «) = -4{V[/ 2 ft ci V • l(r, u) - 6f^ cL Vn (r) • - V x [f h xcT n,V x 3 1 ! ' " ! 



or " n [r) " n [r) 



n Q (r) n (r) n {r) 

+ 3S VM1^1 + 2 Mm : V^VM|^ _ v /iT V • (19) 

where f xcL and f xcT are functions of the local density n (r) and the frequency uj. 

Our result simplifies considerably in two limiting cases. (1) The density is slowly varying 
on the scale set by the wavelength of the perturbation. In this regime we can neglect the 
terms containing Vn , and we recover the former results of GK [0 for the response to a scalar 
potential, and of Ng [16] for the response to a general vector potential. This result can also 



be simply obtained from eq. fpj]) by the substitution n — ► no(r). (2) The velocity field ji/riQ 
is constant in space. This is the case when the static external potential is parabolic and 
the perturbing electric field is uniform, the regime of Kohn's theorem [0 and the HPT || . 
Then all derivatives of j\/n vanish, and eq. (^) reduces to 



ai* c (r» = — 2 {V[f h xci y ■ UrM ~ 5fL L Vn (f) • J -^A. (20) 

By a gauge tranformation, this longitudinal vector potential can be transformed to the scalar 
xc potential proposed by Dobson || to satisfy the HPT. However, for the general case of 
slowly varying no(f) and ai(r, u>), the full expression flT9|) is required. 

In conclusion, our analysis uniquely specifies a local current density functional theory of 
the linear current (and density) response, which becomes exact in the limit of slowly varying 
unperturbed densities and perturbing potentials. The scale on which the variations must 
be slow is set by the smaller of the wave vectors kp and u/vp. Therefore, this theory is 
applicable to the study of high frequency phenomena, such as electromagnetic absorption, for 
which the adiabatic approximation @ is in general not justified. A complete local current 
density response theory, for spatially slowly varying unperturbed density and perturbing 
field, and all frequencies (0 < u> < oo) remains to be developed. 

We acknowledge support from NSF grants No. DMR-9403908 and DMR-9308011. One 
of us (GV) gratefully acknowledges the hospitality of the Institute of Theoretical Physics 
and of the Physics Department at UCSB, where this work was initiated. We thank A. K. 
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